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Abstract 

The 5G wireless network architecture will bring dense deployments of base stations called small 
cells for both outdoors and indoors traffic. The feasibility of their dense deployments depends on the 
existence of a high data-rate transport network that can provide high-data backhaul from an aggregation 
node where data traffic originates and terminates, to every such small cell. Due to the limited range 
of radio signals in the high frequency bands, multihop wireless connection may need to be established 
between each access node and an aggregation node. In this paper, we present a novel transmission 
scheme for wireless multihop backhaul for 5G networks. The scheme consists of 1) group successive 
relaying that established a relay schedule to efficiently exploit half-duplex relays and 2) an optimized 
quantize-map-and-forward (QMF) coding scheme that improves the performance of QMF and reduces 
the decoding complexity and the delay. We derive an achievable rate region of the proposed scheme 
and attain a closed-form expression in the asymptotic case for several network models of interests. It is 
shown that the proposed scheme provides a significant gain over multihop routing (based on decode-and- 
forward), which is a solution currently proposed for wireless multihop backhaul network. Furthermore, 
the performance gap increases as a network becomes denser. For the proposed scheme, we then develop 
energy-efficient routing that determines groups of participating relays for every hop. To reflect the metric 
used in the routing algorithm, we refer to it as interference-harnessing routing. By turning interference 
into a useful signal, each relay requires a lower transmission power to achieve a desired performance 
compared to other routing schemes. Finally, we present a low-complexity successive decoder, which 
makes it feasible to use the proposed scheme in practice. 

Index Terms 


Multihop relay networks, wireless backhaul, physical-layer network coding 
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Fig. 1. Routes obtained via a routing algorithm in an network with multiple sources and a single destination. A dashed circle 
represents a communication range whose radius can be determined as function of a transmit power. Red-dashed lines indicate 
strong interference. Any signal outside the communication range is relatively small and captured as additive Gaussian noise. 


I. Introduction 

To cope with the exponential inerease in mobile data traffie, cellular networks are undergoing 
a paradigm shift from a well-planned deployment of tower-mounted base stations (BSs) to a 
capaeity-driven (possibly ad hoe) deployment of smaller and lower-power BSs ealled small cells. 
This new arehiteeture will eover both the outdoors and indoors seenarios and typically operate in 
higher frequeney bands (mmW) where more speetrum is available [[T|. While conventional BSs 
are typically connected to a high eapacity point-to-point baekhaul network, the same will not 
be possible for small eells due to their dense deployment at often more adverse locations. For 
such networks, the baekhaul is a major bottleneek that can limit the throughput and deployment 
cost 0 . To overcome this problem, baekhaul in 5G needs to adapt to the small eell arehiteeture 
and bring new, more effieient solutions. Towards that goal, 5G baekhaul is expeeted to be 
more integrated with the access network with respect to both technology and spectrum. In 
particular, wireless self-backhaul Q-Q where the same radio speetrum is used for both aceess 
and transport, is an attraetive solution that can substantially reduce the deployment cost, requiring 
only a single radio at eaeh node and allowing for more flexibility in radio resouree alloeation. 
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Fig. 2. Group successive relaying for wireless backhaul networks when L = 2 and K = 3. In each time slot, relays that 
are transmitting are shown in red circles and relays that are receiving are shown in hlack circles. The sources transmit and the 
destination receives in every time slot. Each path contains the L routes for L sources. 


Another key departure from the traditional baekhaul is that, due to the limited range of radio 
signals at high frequeneies Q, a multihop wireless eonneetion may need to be established 
between eaeh aeeess node (AN) and an aggregation node (AgN|^ where an optieal wired or high 
speed radio link eonneetion is available (see Fig. [^. In wireless multihop baekhaul, an aeeess 
node serves not only its own assigned user terminals (UTs) in the vieinity, but also neighboring 
aeeess nodes as a relaying node, in order to forward their data towards and/or from an AgN. 
Therefore, messages are routed (relayed) over multiple wireless hops to reaeh their destination. 

Due to the broadeast nature of wireless medium, interferenee beeomes a main limiting faetor 
on the network throughput for wireless multihop baekhaul network. In Q-Q, interferenee- 
aware routing was proposed and showed to offer a signifieant throughput gain over shortest-path 
routing. A joint routing and resouree alloeation for wireless baekhaul networks was presented in 
0. Beeause it is assumed that eaeh relay deeodes its desired message by treating other signals 
as noise, interferenee-aware routing aims at avoiding inter-route interferenee. It was shown in 
Q that this approaeh ineurs signifieant limitation on network throughput at high load (i.e., the 
number of sourees is large). This result is expeeted sinee it is nearly impossible to avoid all 
inter-route interferenee at high load. Furthermore, beeause the transmission rate on every route 

'Notice that our work can be immediately applied to the uplink of radio access networks (C-RANs) with a multihop backhaul 
network in (To) where sources, relays, and destination in Fig. [T] correspond to mobile stations, radio units, and control units, 
respectively. 
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is determined by the minimum of all link-capaeities on the route, strong interferenee on any one 
link of the route ean drastieally degrade the end-to-end performanee (see Fig. 1). Instead, in this 
case, a more advanced coding scheme that can efficiently manage strong interference instead of 
simply treating it as noise should be considered. We develop such a scheme in this work. 

Our proposed approach is motivated by the view of the wireless backhaul network as a 
multiple-multicast relay network. In such networks, the communication scheme that achieves the 
best known performance is quantize-map-and-forward (QMF) pT| (and a more general scheme 
of noisy network coding (NNC) p^ , p3||). QMF/NNC achieves the cut-set upper bound (on the 
total network throughput) within a constant gap where the gap grows linearly with the number of 
relay nodes. This gap, however, may not be negligible for the system with multihop transmissions 
such as wireless backhaul network operating at high frequencies. In this paper, unlike [ 11 |-p^, 
we assume that each node operates in half-duplex mode, as is the case in many practical systems 


[14|-[16|, due to the problem of saturating the radio front-end when transmitting and receiving 
in the same frequency band and during the same time slot. This makes the considered problem 
more complicated since we need to determine the transmit/receive mode of each node for a 
given time slot. In fact, finding an optimal transmission strategy for wireless multihop backhaul 
networks with either full-duplex or half-duplex nodes is still an open problem. 

In this paper, we present a novel transmission scheme for a wireless backhaul network 
where L sources transmit their independent messages to the destination via {K + 1) hops. 
For each hop, we form a relay stage that consists of 2L half-duplex relays. We refer to such 
network as a layered network. The layered network structure is formed via a routing algorithm 
that we will optimize. In contrast to routing, in our scheme, data from every source message 
is simultaneously forwarded by L relays at each hop. This is achieved by group successive 
relaying and a coding scheme that we develop and that uses the principle of QMF/NNC that 
we develop. In group successive relaying, 2L relays at each stage are divided into two groups 
of equal size so that while the relays in one group transmit their messages to the next stage, 
the relays in the other group receive signals from the previous stage. The role of each group 
is swapped at each time interval (see Fig. |^. In this way, the sources send L new messages 
to the destination at every time slot. Hence, this relaying scheduling is optimal in the sense of 
maximizing the number of data streams (i.e., sources’ messages). Under the group successive 
relaying, we present an “optimized” QMF which ensures a high rate per each message. In the 
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conventional QMF pT| , each relay first quantizes its received signal at the background noise 
level, then randomly maps it to a Gaussian codeword, and transmits the codeword. The random 
mapping can be constructed by random binning and encoding the bin index with a Gaussian 
codeword [|T7|. While DF relays forward one of source messages to the next stage, QMF relays 
forward a function of all source messages (bin indices) as in network coding to the next stage. 
The destination does not explicitly decode the bin indices forwarded by relays. Instead, the 
source messages are decoded jointly with the bin indices, using the time-expanded network in 
Fig. 1^ While NNC scheme does not require time-expanded graph due to a vector quantization, 
it also performs joint decoding over multiple transmission blocks. This requires a joint decoder 
over the entire network in Fig. requiring a high computational complexity at destination. 
The proposed coding scheme (called “Optimized” QMF) is different from the conventional 


QMF [111 in the four-fold: 1) Destination performs forward decoding in which it can start 
decoding a source message after one time slot, using a subnetwork (denoted by Si in Fig. 
and the joint decoding is separately performed for each stage within a subnetwork, whereas it 
is done at once over the entire (time-expanded) network in pTf , [ [T^ . Thus, the proposed QMF 
significantly reduces the decoding complexity and the delay; 2) Destination explicitly decodes 
all relays’ messages (bin indices) and exploits them as side-information in the next time slot, 
which can completely eliminate interference among subnetworks; 3) A quantization level at 
each relay is optimized instead of choosing a background noise level; 4) Vector quantization 
is used as in 0, (H whereas QMF uses scalar quantization. We show that using optimal 
quantization outperforms the other quantization methods that use background noise-level [ |TT| , 
stage-depth [19|,^or classical Wyner-Ziv quantization [20|. Furthermore, both using the stage- 
depth and optimal quantization attain better rate-scaling that the other two methods, i.e., the 
performance degrades logarithmically with K, thereby attaining a larger performance gain as 
K grows. This result shows that optimizing quantization levels at relays substantially improves 
the rate-scaling and is indeed required for multihop transmissions. Furthermore, we show that 
the proposed scheme provides a much higher rate than interference-aware routing based on 
decode-and-forward (in short, MR). The performance gap again increases as the network becomes 
denser (or load becomes higher). We further show that the proposed scheme has a higher energy 


^In stage-depth quantization, the quantization level is chosen at a resolution decreasing with the number of stages K. 
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efficiency compared to MR. The improved energy efficiency comes from enabling each relay to 
collect the signals broadcasted by transmitters (instead of treating them as noise). Based on these 
results, we propose interference-harnessing routing in which the routing criterion is in contrast 
to the routing criterion of interference-aware routing Q-0: the routing metric is chosen to 
“harness” interference, instead to avoid it. Finally, to overcome the drawback of NNC which 
depends on prohibitively complex joint decoding at the destination, we develop a low-complexity 
stage-by-stage successive MIMO decoding, which can be implemented by using a conventional 
MIMO decoding (e.g., zero-forcing, linear MMSE, and integer-forcing receivers) instead of using 
complicated joint decoding. 

The paper outline is as follows. In Section |I^ we describe the proposed scheme and derive 


its achievable rate region. In Section III asymptotic analysis is presented and a closed-form 
rate expression is derived. In Section we develop an interference-harnessing routing suitable 
for the proposed scheme. The low-complexity successive decoder is presented in Section |Vj In 


Section VT numerical results for small networks are provided. Section VII concludes the paper. 


II. The Proposed Transmission Scheme 

We consider a wireless backhaul network where L sources want to send their messages to 
destination with the aid of intermediate half-duplex relays. The proposed scheme consists of 
three parts: i) a routing algorithm performed to establish a {K -f l)-hop 2L-layer network; ii) a 
relay scheduling referred to as group successive relaying applied to assign the transmit/receive 
mode of each half-duplex relay for a given time slot, which is optimal in the sense of maximizing 
the number of messages (simultaneously transmitted to destination); iii) the optimized QMF to 
ensure a high data rate per each message. 


A. Establishing a layered network 

A routing algorithm (e.g., interference-harnessing routing proposed in Section [TV]) is performed 
to establish L routes for L sources. As shown in Fig. the number of hops (or relay stages K) 
of each route is chosen according to the communication range and source-destination distance 
where the communication range can be determined as a function of the transmit power. For 
the convenience of explanation, it is assumed that each route is composed of K relay stages 
(i.e., {K -\- 1) hops). However, the proposed scheme can be naturally extended into the so-called 
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asymmetric layered network where the routes from different sourees ean have different number of 


hops due to the various souree-destination distanees. This ease will be explained in Seetion IV 


To overeome the half-duplex eonstraint, eaeh stage of a route is formed by two half-duplex 
relays, used in transmit and receive modes, such that while one relay transmits its signal to 
the next stage, the other relay receives a signal from the previous stage. Namely, the routing 
algorithm establishes two routes per each source. Therefore, a (K + l)-hop 2L-layer network is 
produced as shown in Fig. [^for L = 2 and K = 3. 


B. Group Successive Relaying 

In group successive relaying, relays are divided into two groups depending on their trans¬ 
mit/receive mode such that 

= {{»?»}£.,■.j = l,...,L} 

Q'l — ■ J — 1. ■ ■ ■ 5 7/} , 

where Ki = [{K — 1)/2J and K 2 = {{K — l)/2], and represents the relay at stage k 
on the 7-th route of source j. At odd (resp. even) time slot (consisting of n channel uses), 
the relays in Qi (resp. Q 2 ) operate in transmit mode and the relays in Q 2 (resp. Qi) operate 
in receive mode. In this way, each source j transmits wj. G {!,... at rate to the 

destination at time slot t, where 7 = 1 for odd time slot t and 7 = 2 for even time slot t, and 
the destination can decode L new messages (wj_^,..., w t-K)- Here, we used two different rates 
Tj 1 and Tj 2 per each source since the odd-indexed and even-indexed messages are conveyed to 
the destination via two disjoint paths: path 1: ((Si,..., S^), (Ri^i,..., R{ : j = 1,..., L), D) 
and path 2: ((Si,..., Sl), (R^^i,..., R^^^ : j = 1,..., L),D). Notice that each path consists of 
L routes for L sources. As shown in Fig. it is assumed that each transmit signal is received 
only at the neighboring relays in the received mode, namely, relays at the same stage on the 
other path and relays at the previous and next stages on the same path. This assumption may 
be reasonable since the strengths of signals from the outside of neighboring relays are relatively 
small thereby being able to be captured by additive Gaussian noise (see Fig. [^. 

Incorporating the group successive relaying into a channel, the discrete memoryless channel 



side-information: 


Destination 


2-stage relays 


1-stage relays 


Sources 


^ \ 


decoded messages: {■^ 3 i^ 4 >w' 2 }i=i 



Fig. 3. Time expanded multihop network with L = 2 and K = 2 where relays in transmit mode are shown at each time slot. 
Red lines stand for the “known” interferences at destination and each arrow actually represents the wireless MIMO channel with 
two inputs and two outputs. 


is described by the transaction probabilities as 

L \K/2] 

/ 7 I 7 1 T 1 T 

F' \yi,2k-l\-^i,2k-2i • • • > ■^i,2k-li ■^i,2ki ■ ■ ■ ) •^i,2k 

j=l k=l 
[K/2\ 

n p 

k=l 

where z = 1 for odd time slot and z = 2 for even time slot, and where and denote 

respective input and output at relay R;’^, and x{ q and denotes source j’s inputs. 

C. Optimized QMF 

We present a novel coding scheme for wireless backhaul networks using group successive 
relaying as the underlying relay scheduling. Accordingly, the channel model in ([^ is considered. 
Fig. 1^ shows the time-expanded network for the case of L = 2 and K = 2. This will be used 
for the decoding of the proposed QMF scheme to deal with cycles in the original network as in 
HD- To clarify our contribution on coding scheme, we compare the proposed scheme (named 
“optimized” QMF) with the previous schemes of QMF, NNC, and SNNC: 


yi,2k 


I i 1 


, X, 


i,2ky 


■^j,2fc+l’ • 


> ^i,2fc+l 


p{yD\xl 


Ki 


(1) 
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i) QMF and NNC in [111, [ [T^ consist of message repetition eneoding (i.e., one long message 
with repetitive eneoding), signal quantization at relay, and simultaneous joint deeoding on the 
reeeived signals from all time slots (say, t = 1,2, ...,T) without explieitly deeoding relays’ 
messages (i.e., quantization indiees). 

ii) SNNC in Q overeomes the long delay of NNC, by transmitting many short messages in 
time slots rather than using one long message with repetitive eneoding. Instead of simultaneous 
joint deeoding, destination performs backward deeoding. After T time slots, the destination starts 
deeoding the souree message and relay’s messages in the slot t, for f = T, T — 1,..., 1, in 
deseending order. Notiee that joint deeoding is separately performed for eaeh slot, instead of 
jointly over all time slots. 

iii) The optimized QMF also uses short messages as in SNNC but performs forward deeoding. 
The time-expanded network is partitioned into T subnetworks, eaeh of them is denoted by St 
for f = 1,... ,T, as shown in Fig. After time slot 1, the destination deeodes the messages 
i.e., souree message and relays’ messages) in the subnetwork 5i, and then, after time slot 2, 
deeodes the messages in ^2 by exploiting the previously deeoded messages as side-information. 
Due to the use of sueeessive relaying, the destination ean eompletely know the interferenee seen 
at relays (depleted by red arrows in Fig. sinee they are fully determined by the previously 
deeoded relays’ messages. In this way, the destination ean deeode the all souree messages 
for f = 1,..., T, in that order. Thus, the forward deeoding has a lower deeoding delay than the 
baekward deeoding in Q. Furthermore, it enables to find an optimal quantization level at eaeh 


relay, whieh will be explained later in Seetion 111 


Taking our deeoding method (forward deeoding) into aeeount, we are able to produee the 
simplified ehannel model illustrated in Fig. where we introdueed the notation of ^ sinee the 
known interferenee ^ is different from x] ^ with respeet to a time index. Sinee this model ean 
be applied to deeoding of any souree message regardless of time indiees, we will use this model 
to derive an aehievable rate of the proposed seheme. For the ease of exposition, we introduee the 
notation • • • ’• • •) ^,k) that denotes the known interferenee (i.e., 

side-information) at the destination. That is, the destination eompletely knows (Xj i,... ,Xj i<') 
when deeoding souree messages eonveyed via path i. 

From now on, we will explain the eneoding/deeoding proeedures foeusing on the ease i = 1 
(i.e., for the network model in Fig. Q. The same proeedures apply to the ease i = 2. The 




side-information 

(11.1.11,2) 
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Fig. 4. Equivalent simplified model of 3-hop 2-layer network where = (^2 ir-n ■ ■ ■ 1 ^2 ir-n m • • • 1 it) 

denotes the known interference at destination. 


encoding/decoding procedures are described as follows. 

Codebook generation: Fix input distributions as 

2 K L 

]l]l]lpKk)pifi,k\yik)- (2) 

i=l k=0 j=l 

Randomly and independently generate 2"'^i codewords xloK) of length n indexed by wl E 
{1,..., 2’^^*} with i.i.d. components ~ p{xIq). Randomly and independently generate 

• • • j 

codewords of length n indexed by E {1,..., 2"^’'^'=} with i.i.d. components ~ 

‘'3 ■ '■3 

Independently generate 2"^^*^ codewords y\j3^) of length n, indexed by z/ G {1 ,..., 2"^^hfc}, 
with i.i.d. components ~ PiVik)- The quantization codewords are randomly and independently 
assigned with uniform probability to bins. We denote the bin by with E 

2 <^}. 

Encoding: At every odd time slot, each source j transmits a new message w{ G 2”^i} 

by sending the codeword x^ Q(yy{). Each relay R{ ^ observes the ^ in receive mode. The 
relay quantizes its observation into a quantization codeword yj ^ and finds a bin-index i\j^E 

j 

{1 ,..., 2’^'^hfe} such that the corresponding bin contains the quantization codeword yifc- In 
transmit mode, the relay transmits the bin-index by sending the codeword ^). 

Stage-by-stage successive decoding: At every time slot, the destination proceeds to decode 
{^\ ^,... i^ik) (i-e., the messages of relays at stage k) for k = K,... ,1, in that order, and 
then decode source messages (w|,... ,w^). Notice that the destination knows (Xi i,... 
as side-information since as explained before, they are completely determined by the decoded 
messages in the previous time slot. 
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From the received observation the destination can reliably decode k, ■ ■ ■ 

e 7^ (p( 2 /d|x}_^,. . , (cx),..., cx))) , 

where the above rate region is equivalent to the MIMO MAC capacity regioij^and where 7^(-) 
is given in Definition below. Using the decoded bin-indices the side- 

information Ii^K, the destination performs joint typical decoding to decode • ■ ■; 

namely, it finds unique message {^\ K-l^ ■ ■ ■ , ^yK-l) such that 

for some ^ W{,k)^ 

where Te"'\Xi^o, Yi, denotes the jointly typic3.1 set for fixed 

..., rfx), the channel for decoding the above messages can be seen as the multiple access relay 
channel (MARC) consisting of L sources, L intermediate relays, and one destination. Specifically, 
L sources transmit their signals to L relays via wireless channel defined as 

L 

i=i 

and L relays send their messages to the destination via noiseless links of capacities (rtK, 

From Lemma [T] and Definition [T] below, we have: 

• • • y^yK-i) ^ {yi,K\^i,K-iy ■ ■ ■ y ^i,K-iy^i,K) , 

With the same argument, the destination can decode {j\^K- 2 T ■ ■ ^ ^i,k- 2 ) using the decoded 
bin-indices {i] j^_i,..., if j^_i) and the side-information Ii^k-i if 

• • • 5 ^ {yi,K-l\^l,K-2y • • • 5 ^l,K- 2 ^^l,K-l) ) • • • 5 ‘^yK-l) 

In general, the following rate-constraints should be satisfied to reliably decode all relays’ mes¬ 
sages: For k = K,... ,1, 


XyK)- 


^In the MARC, if the link-capacity between each relay to destination is infinity, then the destination can observe the received 
signals at the relays without any distortion. Thus, this model is equivalent to the MIMO MAC channel. 
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Notice that the above rate regions ean be derived using K + 1 recursively defined ehannels. 
Finally, the destination ean reliably deeode L new messages (w},..., wf) using the deeoded 
bin-indiees ..., and the side-information Ii i if 

The exaetly same analysis ean be applied to the ease of z = 2. From the above analysis and 
letting r- = r^g for the eonvenienee of notation, we obtain: 

Theorem 1: For the wireless baekhaul network in Fig. a rate tuple (rj,..., rf, ..., ) 

is aehievable if 

L 


e 7^ I Wp{y\ 

\j=i 




) • • • 5 '^j,A:+l) I ! 


fox k = K^K — 1,0, with initial value = cx), j = 1 ,..., L, for any input distributions 

defined in ■ 

Lemma 1: Consider the MARC defined as 

Ylp{yik+iHk, ■ ■ , {rlk+i, ■ ■ • ,C>+i) 

Kj=l 

When Xjfc+i is known to the reeeiver, a rate tuple is aehievable if for all C C 


^Ik < min [<fc+i - ^ {"^lk+i'^%+i\^lk^ • ■ ■ > 
iec j£S 

+ I {(Xl, : e £); : j e S‘)\(Xl, : j e . 

for any input distributions that faetors into Ylj=iPi^ik) Y[j=iPiyik\yik)- 

Proof: The proof almost follows p8l Appendix 1]. The eonsidered model reduees to the 


model in [18| by setting Xj ^+i = f. Sinee the known interferenee Xj ^+i ean be eompletely 
eaneeled at the destination, it does not ehange the aehievable rate region. This eompletes the 
proof. ■ 

Definition 1: We denote the rate region of Lemma by 

^ {\\p {yik+iWi,k^ • • • > ^ik^^k+i) , (rjfc+i,..., 


(3) 
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Remark 1: A classical compress-and-forward (CF) p0[ can be considered as a special case 
of QMF in the sense that a distortion (or quantization) level at each relay is chosen according 


to Wyner-Ziv distortion [20|. That is, the quantization level at relay is chosen such that 


= ri,,. 


(4) 


Then, the destination can avoid the complexity of joint typical decoding and just use classical 
successive decoding. Namely, the destination can find quantized observations using the decoded 
bin-indices and side-information. Then, it just performs “classical” MIMO decoding based on 
quantized observations, which significantly reduces the decoding complexity at the destination. 
Thus, using Wyner-ziv quantization may be better choice in practice as long as it provides a 
satisfactory performance. 

Remark 2: In Theorem we derived an achievable rate region of the proposed scheme. 
However, it is not so clear how to optimize relays’ rates subject to their achievable rate region 
for maximizing a sum rate or a symmetric rate. This is non-trivial optimization problem and 
is interesting subject of a future work. In the next section, we partially solve this problem by 
restricting our attention to symmetric Gaussian networks. 


III. Gaussian Wireless Multihop Backhaul Networks: 

Asymptotic Analysis and Closed-Form Rate Expression 

We consider a Gaussian wireless backhaul network where destination is equipped with L 
receiver antennas and accordingly, the L-dimensional observations are denoted by • • •, 


). Let Hjfc denote the Lx L channel matrix with inputs (xj^,..., and outputs 

for k = 0,1,..., iT. Let denote the L x L intra-path interference 


Ai,K+i 




and outputs (yj^,..., yJ^J for /c = 1 ,..., AT - 1. 


channel matrix with inputs 
Also, let denote the Lx L inter-path interference channel matrix with inputs • • • ^xf^) 
and outputs (y^ ,,■■■, y^.) for fc = 1,..., K. Lrom Theorem 1 we can see that the performance 

- 7/ f fC - X ^ rC 

of the proposed scheme is independent from the intra-path and inter-path interference matrices 
Gj fc and Sj^^. Lor the Gaussian channel. Lemma is simplified as 


Lemma 2: Lor the Gaussian MARC defined by (Hj ^, 
is achievable for all £ C {0 ,..., L — 1}, 


,G,fc+D • • •) G,fc-i-i 


)), a rate tuple 


i.k") ' 


£ec 


5C[1:L] ^ 
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for some quantization level where represents the ehannel sub-matrix to eontain 

the rows of with their indiees belong to 5 C [1 : L], ^ represents a diagonal matrix 

whose j-th diagonal element is Ptx/(1 + k)' denotes a transmit power at eaeh node.B 

We evaluate the performanee of the proposed seheme for two extreme seenarios of network 
strueture ealled sparse and dense networks. They are elassified aeeording to the strueture of 
In a sparse network, the number of interfering nodes is eonstant as the number of nodes L grows 
(i.e., Hj fc forms a band matrix) while in dense network, it inereases linearly with L (i.e., Hj ^ 
forms a “full” matrix). Although we limit our attention to Gaussian networks, the eomputation of 
the aehievable rate in Lemma is generally diffieult sinee it involves a eomplieated eombinatorial 
optimization. In order to make a problem manageable, we eonsider the partieular struetures of 
Hj fc that eaptures the features of sparse and dense networks, respeetively. They are deseribed 


in Seetion III-A and III-B, respeetively. Further, we let L —)■ cx) (i.e., asymptotie analysis) and 
foeus on the aehievable symmetrie rate of r. 


A. Sparse Networks 

In this seetion we assume that eaeh transmit signal is reeeived at its desired relay and two 
neighboring relays, i.e., has the form of tri-diagonal strueture with (a, 1, a), where a denotes 
the relative interferenee strength. Namely, other interferenees are relatively weak and are eaptured 
by additive Gaussian noises. This model is well-known as Wyner model p2| . Obviously, this 
model is the form of sparse network in that the number of interfering nodes is fixed by 2, 
independently from the number of users L. The transmit power at eaeh node is assumed to 
be Ptx = SNR. Hereafter, we restriet ourselves to ehoose relays’ rates as for all i 

and j, whieh is reasonable due to the symmetrie strueture of Hj ^. Then, by symmetry and 
eoneavity, this limits the sum-rate inequality to be dominant in Lemma and by ehoosing the 
same quantization levels at eaeh relay (i.e., Qk = Qlk i)’ we have: 

= I .SrA . 

( 6 ) 


A remarkable result of [18| is that in the limit of L — >■ cxd, Q ean be simplified to Q below. 
Consequently, we get the following: 
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Corollary 1: The proposed scheme achieves the symmetric rate (per user) r = tq to satisfy: 


rk-i = F{x*), 


(7) 


for k = K + 1,... ,1, where 

F{x)= [ log(l + SNR(l-2-^)(l + 27cos(27r0))2))ci0, 

Jo 

and X* is the solution of the equation F{x) = Vk — x with the initial value tk+i = oo. 


For the comparison, we consider the MR where each relay decodes its desired message 
(forwarded by its own route) by treating all other signals as noise. In order to derive its achievable 
rate, we assume that Gi^k = 0 and 8*,^ = 0. Thus, the achievable rate of the MR is optimistic 
since in practice, they may not be zero matrices. Under this assumption, the achievable symmetric 
rate of the MR is given by 


= log 1 + 


SNR \ 


( 8 ) 


l + 2a2SNRy ■ 

Fig. 1^ shows the achievable symmetric rate of the proposed scheme for various quantization 
levels at relays. We first observe that using optimal quantization outperforms the other cases 
as background noise-level (i.e., Qk = 1) GD’ stage-depth (i.e., Qk = K) p9| , and Wyner-Ziv 
quantization (obtained from Q) p0| . As anticipated, the optimal and stage-depth quantization 
provide a larger performance gain due to the improved rate-scaling. When K > A, the noise- 
level and Wyner-Ziv quantization show worse performance than the MR and hence optimizing 
a quantization level plays an important role for multihop networks. 


B. Dense Networks 

In this section we assume that Hj ^ is an i.i.d. matrix whose entry has a complex Gaussian 
distribution with zero mean and unit variance. Note that the transmit power at each node is 
assumed to be Ptx = SNR/L. If we choose the Ptx = SNR as in the sparse network, the received 


power at each relay goes to infinity, thereby yielding an infinite rate. As done in Section III-A 
we limit our attention to the symmetric choice of for all j and i. With the same 

argument in Section III-A, the sum-rate inequality in Lemmais dominant and is represented 
as in ([^ with Ptx = SNR/L. A closed-form expression of (|^ has been derived in [23| by 


using the fact that as L —)■ cx), the problem becomes symmetric although the channel is “full” 
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10 


Number of stages (K) 


Fig. 5. SNR = 20dB and INR = 15dB (e.g., a = -y/lNR/SNR = 0.56). Performance of the proposed scheme for various 
quantization levels as a function of the number of stages K. 

and non tri-diagonal as in the Wyner model. The elosed-form expression was obtained from the 
asymptotic random matrix theory and the submodular structure of the rate expression, which 
equal to Q below. Consequently, we get the following: 

Corollary 2: The proposed scheme achieves the symmetric rate (per user) r = tq to satisfy: 



(9) 


for fc = + 1,..., 1, with initial value Qk+i = 0 and tk+i = C) 0 , where 



( 10 ) 


Since the achievable rate in Q is the minimum of two terms, where the first is an increasing 
function of Qk and the second is a decreasing function of Qk, the optimal value of is attained 
by solving 


rfc - log(l + 1/Qfc) = C (SNR/(1 + Qfc)) • 


f{Qk) =n- log (1 + 1/Qk) - C (SNR/(1 + Q,)), 


Letting 
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Fig. 6. LSNR = 20 dB. Performance of the proposed scheme for various quantization levels as a function of the number of 
stages K. Notice that the performance of the MR goes to zero as L —> oo due to the impact of severe interference. 

we can find Qk,mm = 1/(2'’'' - 1) and Qfc.max = (1 + SNR)/(2^'= - 1) such that f{Qk,mm) < 0 
and /(Qfc,max) > 0. This is because QA;,min makes the first term of the minimum in (|^ zero 
and Qfc,max is the Wyner-Ziv quantization, which makes the second term to attain the minimum 
in Q. Using bisection method, we can quickly find an optimal quantization level QA:,opt- This 
value will be used to plot the performance of the optimized QMF. 

Fig. 1^ shows the achievable symmetric rate of the proposed scheme for various quantization 
levels. The performance trend is similar as in the case of sparse network. The most remarkable 
observation is that the proposed scheme can achieve almost the same performance as in Fig. 
with a lower transmit power (per each node) as L —)■ cx). 

IV. Interference-Harnessing routing 

The proposed scheme can be applied together with any multihop routing algorithm (e.g., 
interference-aware routing Q). In this section, we provide a routing criterion that further im¬ 
proves the performance of the proposed scheme. From Figs.|^and[^ we observe that the proposed 
scheme can achieve almost the same performance in both dense and sparse networks. In dense 
network, however, the transmit power consumption is reduced proportionally to the number of 
users L (i.e., Ptx = SNR/L). In other words, the proposed scheme has a higher energy efficiency 
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Fixed route (established at the previous iteration) 






Iterative Algorithm: 

fixed the new route and update the previous route 
(d) step 4 


Fig. 7. Illustration of an interference-harnessing routing to establish the two routes in path i where “candidate relays” denotes 
the relays within the communication range from a transmitter. 


achieving a desired rate performanee with a lower transmission power as the network beeomes 
denser. The improved energy effieiency eomes from enabling eaeh relay (by using QMF) to 
collect the signals broadeasted by transmitters. Any interfering signal that is reeeived at the 
relay will be forwarded through QMF and treated as a useful signal at the destination. For 
this reason, QMF performs better when the network is dense and interferenee is stronger. This 
motivates us to propose interference-harnessing routing in whieh the routing eriterion eontrasts 
the routing eriterion used in the interference-aware routing. Interferenee-harnessing routing refers 
to a family of routing algorithms in which the metrie is ehosen to exploit interferenee. This can 
be done by, for example, ehoosing the metrie that maximizes an aehievable rate between every 
two conseeutive relay stages. For any ehoice of relays at the two stages, this rate corresponds 
to a rate in a MIMO ehannel (see Seetion |V]) and ean thus be ealeulated, as we will specify in 
this section. An approximation of such metric would simply choose relays at eaeh stage sueh 
that the signal power received from the previous stage at each of these relays is maximized (see 
Remark [^. The performance gain of interference-harnessing algorithm eompared to the MR is 
demonstrated via simulation result in Seetion 

For the interferenee-harnessing routing, one efficient algorithm ean be developed using the 
iterative algorithm in Q by properly modifying the routing eriterion. In Q, the algorithm 
establishes one route at a time while keeping other previously established routes fixed, and 
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repeats the proeess until negligible improvements in the sum throughput ean be made. The 
routing eriterion maximizes eaeh link-eapaeity on the route whieh is eomputed by taking into 
aeeount interferenee from all other routes. For the interferenee-hamessing routing, on the other 
hand, we seleet a relay (at stage k) to maximize the MIMO eapaeity defined by two eonseeutive 
stages k — 1 and k. We provide a deseription of this algorithm for the network with two sourees 
using Fig. Generalization to the ease with more sourees is straightforward. We use the notation 
Cmimo({Ri, R 2 }, {R 3 ,R 4 }) to denote the MIMO eapaeity indueed by two transmitters {Ri,R 2 } 
and two reeeivers {R 3 ,R 4 }. Without loss of generality, we foeus on the path 1 (i.e., i = 1 in 
Fig. 0. The algorithm proeeeds as follows: 1) Establish the route from Si to D so that the 
number of hops is minimized and eaeh link-eapaeity along the route is maximized (namely, the 
reeeived power is maximized). This proeess elosely follows the step that establishes the initial 
route in 2) For the fixed route from Si to D, establish a new route from S 2 to D, as illustrated 
by the steps in Fig. |7] and deseribed below: 

• Find the relays within a eommunieation range from the S 2 that have not been already ehosen 
for other routes. The eommunieation range is a design parameter to be determined as a 
funetion of a transmit power where the transmit power should be larger than a threshold to 


satisfy network eonneetivity [211. Let 71 = {Ri,i,..., Ri,|ri|} denote the set of sueh relays 
(i.e., “eandidate relays” of the first stage). For the j-th relay in 7), we ean eompute the 
MIMO eapaeity (i.e., routing metrie) Cmimo({Si, S 2 }, {R} 1 , Riy}) where R}^ denotes the 
first relay on the route from Si to D. With Si, S 2 , and R} ^ fixed, we ehoose a relay Rij 
to maximize the MIMO eapaeity: 


j* = argmax^.g[i^|7-^|]CMiMo({Si,S2},{Ri_i,Rij}). 

• Let R^ ^ = Riyj. Find the eandidate relays of the seeond stage, denoted by 72 = {R 2 ,i, • • •, 
R 2 ,|r 2 |}- Then, ehoose a relay R 2 J to maximize the routing metrie: 

J 2 = argmax^.gji^|7-^|jCMiMo({R},i,R?,i},{Ri,2,R2,j}), 

where R} 2 denotes the seeond relay on the route from Si to D. 

By setting R^ 2 = R 2 J 2 ’ establish the route from S 2 to the D. Then, for the fixed route 

from S 2 to D, update the route from Si to D. Repeat this proeess until negligible improvements 
in the sum throughput ean be made. 
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(a) Interference-haraessing routing 


(b) Interference-aware routing 


Fig. 8. The opposite trend of the routes selected by interference-harnessing and interference-aware routing. Blue-colored relays 
establish the path 1 and red-colored relays establish the path 2, respectively. 


Notice that since the performance of the proposed scheme is independent of the inter-path 
interference as explained in Section II-C[ we first apply the above procedure to find a path 1 
and then apply the procedure to find a path 2 using only the relays not included in the path 1. 

Remark 3: (Energy harvesting) In the example of Fig. let Hj denote the 2x2 channel matrix 
between the two transmitters {Si, S 2 } and the two receivers (RJ 1 , Riy} for some Rij G 7i, i.e., 


H 


j 


hii hi2 
hji hj2 


where hu denotes the channel from Si to R| 1 along the established path 1 and hji denotes the 
channel from Sj to a candidate relay Rij G 71. Assuming that the transmit power is SNR, we 
have: 


Cmimo ({Si, S2}, {R},i, Riy}) = logdet (I + SNRH.H^) 

< log(l + + |/fi2p)SNR) + log(l + {\hji\‘^ + |/fj2p)SNR), 

where (a) follows from the Hadamard’s inequality and equality is achieved if the two columns 
of Hj are orthogonal. Thus, the obtained upper bound is maximized by choosing a relay Ri ^ to 
maximize SNRTj = + |/ij 2 p)SNR. We observe that SNRr^ is proportional to the power 

of the signal received from Si and S 2 at the Rij. We can use SNRrj as a routing metric 
approximating the MIMO capacity routing metric in the above iterative algorithm. This also 
implies that the interference-harnessing routing tends to choose a relay such that the signal 
power received from the previous stage is maximized. 
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Remark 4: (Minimizing end-to-end delay): For the interference-aware routing (based on DF), 
confining to nearest-neighbor transmissions maximizes the network throughput by mitigating the 


impact of inter-route interference [24|. However, this approach increases the number of hops to 
reach destination thereby yielding long end-to-end delay. For the interference-harnessing routing, 
on the other hand, longer one-hop transmission (i.e., using a higher transmit power subject to a 
transmit power constraint) increases the network throughput by decreasing the number of relay 
stages K. This is due to the fact that, when relays use the optimized QMF, the throughput 
degrades as K grows (see Figs. [^and[^. Thus, the interference-harnessing routing can be also 
more suitable for the systems with a delay constraint due to shorter end-to-end delay compared 
to the interference-aware routing. 

Remark 5: (Potential Gain of the interference-harnessing routing in multiple AgNs): A natural 
extension of our work is to consider wireless backhaul networks with multiple AgNs (say, M 
AgNs). We can define a subnetwork i consisting of AgN i and the associated sources and relays 
for z = 1,..., M. Each subnetwork can be established via interference-harnessing routing. Then, 
the proposed transmission scheme can be applied to each subnetwork separately. In this case, 
there exist inevitable interference caused by the relays associated with different subnetworks. 
We refer to such interference as inter-network interference. Due to the use of the interference- 
harnessing routing, each subnetwork spans a narrow area over the entire network since in order 
to exploit interference, the routes tend to be chosen as closely as possible (see Fig. (a)). 
For the interference-aware routing, on the other hand, each subnetwork spans a wide area over 
the whole network in order to avoid inter-route interference (see Fig. (b)). Therefore, the 
interference-harnessing routing is much more efficient in avoiding the inter-network interference 
than the interference-aware routing. We expect that, when applying the proposed scheme to 
multiple AgNs, interference-harnessing routing further improves the performance compared with 
interference-aware routing. Moreover, we can employ the interference-harnessing routing to 
establish each subnetwork and the interference-aware routing to avoid inter-network interference. 
The algorithm establishes one subnetwork at a time while keeping other previously established 
subnetworks fixed, and repeats the process until negligible improvements in the sum throughput 
can be made. For the fixed subnetworks z for z G {1,..., we can establish a subnetwork 

j as follows: 

• Perform the interference-aware routing to establish the first route of a subnetwork j, where 
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each link-capacity on the route is computed by taking into account interference from all 
other subnetworks. This process can avoid inter-network interference. 

• Given the first route, perform the interference-harnessing routing to establish the subnetwork 

j- 

Remark 6: (asymmetric layered network): One may be concerned that routes that roots from 
different sources will have different number of hops due to the various source-destination 
distances. We refer to such a network as an asymmetric layered network. This issue can be 
addressed by grouping the sources that are closely located and serving them simultaneously, 
which can be viewed as user scheduling. Furthermore, such user scheduling can maximize 
the “interference-harnessing” gain since it is likely to produce a path such that the relays in 
each stage are closely located and hence each relay (using QMF) can collect more broadcasted 
energy. We emphasize that although routes have different number of hops, the proposed scheme 
can be applied to such asymmetric layered networks as follows. Consider the example of Fig. 
where the source i communicates with the destination with (Ki + 1) hops. For the asymmetric 
layered network, the destination performs the stage-by-stage successive decoding in Section 
in order to decode relays’ messages as well as sources’ messages. Until decoding the messages 
of relays at stage 3, the destination follows the same procedures as in Section |I^ Then, using the 
£^ 3 , f'fs) and the side-information Jj^ 3 , the destination can decode the two relays’ messages 
( 4 ^ 2 ) 2 ) the source 2’s message W 2 . Similarly, using the ( 4 \ 2 )^f, 2 ) the side-information 

Xj 2 , the destination can decode the relay’s message 4,1 and the source 3’s message W 3 . Finally, 
the destination can decode the source I’s message yy^ using the and the side-information 
Xj 1 . This example shows that the proposed scheme can be naturally applied to the asymmetric 
layered networks. 


V. Practical Construction of the Proposed Coding Scheme 


It was shown in Section III that the proposed scheme outperforms the MR. Recall that, in the 
proposed scheme, destination must perform joint decoding to decode relays’ messages as well 
as sources’ messages. It is noticeable that the joint decoding here is separately performed for 
each stage while it is done at once over entire messages in GD. p^ . Nevertheless, the joint 
decoding typically requires higher complexity and makes it harder to design a practical coding 
scheme for the proposed scheme. Although some progresses of practical joint decoding have 
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Fig. 9. Equivalent simplified model of the asymmetric layered network with Ki = 4, K 2 = 2, and K 3 = 3, where Ki 
represents the number of relay stages so that the source i can reach the destination. 


been made in [25|, [261 for small one relay network, we need further investigation to extend the 
previous works into the considered multihop network. 

In this section, we develop a low-complexity decoder named stage-by-stage successive MIMO 
decoder where a conventional MIMO decoding is applied to each stage instead of joint decoding. 
First of all, to avoid the joint decoder at the destination, we resort to the Wyner-Ziv quantization 


in which each relay chooses the quantization level such as 


1 + Etil^llPSNR 


QU = 


2\k 


( 11 ) 


where hfl denotes the (j,£)-th element of and denotes the transmission rate of relay 

1 this particular choice of quantization level ensures that the 
destination can find a unique quantization sequence from the decoded relay’s message and 
side-information. 

Using the above fact, the destination performs the stage-by-stage successive MIMO decoding, 
i.e., the following procedures are repeated for k = K 1, K1, in that order: 

MIMO decoding at the stage k: The destination has the side-information Xi^k and 
(i.e., decoded relays’ messages at stage k -\- 1). 

• Using the side-information {&ik}j=i and Xj^, the destination can find unique quantized 
sequences {y^ Jy=i: 


As noticed in Remark 



+ Y 9fUik+i + Y +^i,k+^i,k^ 


e=i 


e=i 

\ _ 


i=i 
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and where z^^i, denotes an i.i.d. Gaussian random variable ~ 

• After removing the known interferenee Xj fc, the destination ean deeode relays’ messages 

the resulting “quantized MIMO MAC” ehannel, given by 


- 1 




+ 


^2,/c 



9L \ 


^i,k + 


where — li k- Equivalently, we have the matrix form as 

- 2 ; K - 2 ) K ^ 

^j }- = “f (12) 


where eaeh row of a matrix eonsists of n-dimensional veetor. 


Sinee eaeh stage k eonsists of the quantized MIMO MAC defined in ( [12] ), many approaehes 
to MIMO deeoding ean be applied. Clearly, the best performanee ean be attained by maximum 
likelihood (ML) deeoding, whieh aetually aehieves the information-theoretieal rates of using 
Wyner-Ziv quantization in Figs.j^andj^ The eomplexity of this approaeh, however, is exponential 
in the produet of the eoding bloeklength n and the number of reeeiver antennas (i.e., L in our 
ease). The eomplexity of ML deeoding ean be signifieantly redueed through the use of sphere 
deeoding algorithms 


|, [ |28| |. Rather than proeessing all the observed signals from the antennas 
jointly, one simple and widely-used approaeh is to separate out the transmitted data streams 
using linear equalization and then deeode eaeh data stream individually sueh as zero-foreing 
(ZF), linear minimum mean-squared error (MMSE), and integer-foreing (IF) reeeivers [ |29| . The 
ZF reeeiver (a.k.a., deeorrelator) inverts the ehannel matrix so that eaeh data stream ean be 
reeovered via a single-user deeoder. The MMSE reeeiver performs the same operation exeept 
with a regularized ehannel inverse that aeeounts for possible noise amplifieation. Both of these 
arehiteetures permit the use of powerful point-to-point ehannel eodes (e.g.. Turbo eode, LDPC 
eode, and Polar eode) that ean aehieve high data rates at praetieally-relevant SNRs. Reeently, IF 


reeeiver was presented in [29| where the deeoder first eliminates the noise by deeoding a linear 
eombination of interfering data streams and then eliminates interferenee between data streams 
in the digital domain. This seheme is based on the faet that eaeh data stream is drawn from 
the same lattiee eodebook whieh ensures that any integer eombination of eodewords is itself a 
eodeword, and thus deeodable at high rates. Further, low-eomplexity eoding frameworks based 


on QAM modulation and non-binary linear eodes haven been proposed in [30|, [32|. 
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In the next, we attain the performanee of the low-complexity scheme with linear receivers. In 
this case, the receiver (i.e., destination) applies the equalization matrix e to obtain 




(13) 


where Bj ^ will be specifically defined according to linear schemes. From this, we can derive 


an achievable rate ^ as a function of B 


'i,k‘ 


^ik = 


(l 

SNR 


) 

V ' Eti l&gl 

n^ + Qik 

) + SNR 1 

ui i,e 1 

2 1 


where denotes the j-th column of and ..., ) denotes the j-th row of 

Bj^fc. With initial value QIk+i = 0 ^ ^e can derive all relays’ rates at stage k for 


k = K,... and sources’ rates recursively. From |29|, we classify the equalization matrix Bj 


i.k 


H,; .RH ^-1 


according to linear schemes: i) Bj^^ = Hjfc foi" ZR receiver; ii) Bj^^ = ^ -i- n.j^fcn.• 

for linear MMSE receiver, where ^ denotes a diagonal matrix with SNR/(1 -f as its 


j-th entry; iii) Bj^fc = ^ for IF receiver, where e is an 

integer full-rank matrix and is optimized as a function of channel matrix, and Qj ^ denotes a 
diagonal matrix with its j-th diagonal element SNR/(1-|-Q^^). If we choose A = I, this scheme 
reduces to linear MMSE receiver. The key step underlying this approach is the selection of an 
integer matrix A to approximate the channel matrix Hj ^. Although finding the optimal Aj ^ 
has a worst-case complexity that is exponential in L, this search only needs to be performed 


once per coherence interval. In practice, efficient approximation algorithms (see [311 and 
for details) can be used to find near-optimal Aj ^ in polynomial time. 


VI. Numerical Results: Non-Asymptotic Case 

In this section, we evaluate the performance of the proposed low-complexity scheme in 
Section |V] for non-asymptotic cases (i.e., small number of users). From this, we can identify 
the practical feasibility of the proposed scheme. We first investigate the performance of the 
proposed scheme for various receiver architectures namely ML, ZF, MMSE, and IE receivers. In 
our simulation, we consider the dense network with i.i.d. channel matrix H* ^ where each entry 
has a complex Gaussian distribution with zero mean and unit variance. Notice that in Eig. [T^ 
ML receiver achieves the information-theoretical achievable rate. Among the linear receivers, 
IE receiver shows the best performance achieving about 1-bit gain over MMSE receiver. Eor 
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Number of stages (K) 

Fig. 10. SNR = 30dB and L — 4. Performance comparison of the proposed low-complexity scheme for various receiver 
architectures as a function of the number of stages K. 



Fig. 11. L = 4 and K = ?>. Performance comparisons of the proposed scheme and multihop routing where the low-complexity 
decoding (with ML) is used for the proposed scheme. 


< 3, IF receiver can provide a satisfactory performance within 1-bit from the ML performance, 
significantly reducing the decoding complexity over ML receiver. Thus, IF receive can be a good 
candidate for a practical system where the number of stages is limited due to a delay constraint. 
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Fig. 12. SNR = 20dB. Performance comparisons of the proposed scheme and multihop routing where the low-complexity 
decoding (with ML) is used for the proposed scheme. 


For K > 4, however, a low-eomplexity ML reeeiver (e.g., sphere deeoding) should be eonsidered. 

We next evaluate the performance of interference-harnessing algorithm and compare it to the 
performance of interference-aware routing. In our simulation, we consider the multihop wireless 
network depicted in Fig. with L = 4 sources and K = 3 stages of relays. We define the cluster 
(i.e., a square in Fig. whose area depends on the transmit power. It is assumed that each relay 
can only communicate with the relays within eight neighboring clusters and its own cluster 
(i.e., communication range). Any interference outside the communication range is captured by 
additive Gaussian noise. We further assume the symmetric channel model for which a channel 
coefficient (within the communication range) is defined as SNRexp(j6') where SNR captures 
the distance-dependent path-loss and transmit power, and 9 ~ Unif[0, 27r) denotes a random 
i.i.d. phase. In our simulation, we assume that the whole network is divided into Ax K disjoint 
clusters as shown in Fig. [^with K = 3, where K captures the source-destination distance (i.e., 
determine the network area). It is also assumed that each cluster contains the Uc half-duplex 
relays with Uc > L. Notice that for a given K (i.e., network area), the network is denser and 
denser as L grows, since the network includes more nodes. Hence, the parameters K and L 
control the network area and the network density, respectively. 

For such network, the interference-harnessing and interference-aware routing are used to 
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establishes the L routes (per eaeh path) for the proposed seheme and the MR, respeetively. 
That is, the iterative algorithm in Q is applied to both oases with different routing metrios. 
For the interferenee-harnessing routing, we employ the simple routing metrie in Remark in 
whieh we only need to oompute the reoeived power at eaeh oandidate relay and ohoose a relay 
with the maximum reeeived power. In the example of Fig. the relays ehosen by the algorithm 
are shown in the blue and red eireles. This oan be easily extended into a general L due to the 
assumption of Uf. > L. To be speoifio, the interferenee-harnessing routing ehooses the 6 olusters 
as in Fig. (a) where eaeh aotive oluster oontains the L seleoted relays while the interferenoe- 
aware routing ehooses the 12 elusters as in Fig. [^(b) where eaeh aetive oluster eontains the L/2 
seleoted relays. For the MR, there are two types of relays aeoording to the number of reoeived 
interferenoe signals where the relays in the top or bottom route reeeive the (|L — 2) interferenoe 
signals and the relays in the middle routes reeeive the (2L — 2) interferenoe signals. Therefore, 
the aehievable symmetrio rate of the MR is determined by the minimum of all message rates, 
given by 

tmk = log (1 + SNR/(1 + (2L - 2)SNR)). 


Fig. shows that the proposed seheme outperforms the MR having a larger performanee gain 
as SNR inoreases. This is beoause the strong interferenoe limits the performanee of the MR 
while it further improves the performanee of the proposed seheme. From Fig. we observe 
that in eontrast to the MR, the performanee of the proposed seheme is improved as the number 


of users L inoreases, whieh is well-matohed to the asymptotio result in Seotion III-B 


VII. Conclusion 

We presented a novel transmission seheme for wireless baokhaul networks. In our seheme, L 
sourees transmit their messages over multiple hops using L layers at eaeh hop. In eontrast to 
routing, in our seheme, message from every souroe is simultaneously forwarded by all L relays 
at every hop. This is aohieved by developing QMF/NNC by whieh a relay forwards a funotion 
of all souree messages simultaneously. 

We optimized both the relay seleotion as well as the QMF seheme developed by the relays. 
In partieular, we proposed the optimized QMF that performs optimal quantization of reoeived 
signals at the relays. As we showed, the optimized QMF signifioantly outperforms other forms 
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of QMF as well as interference-aware routing (based on decode-and-forward) in both sparse and 
dense networks. Furthermore, the obtained performance gain increases as the network becomes 
denser (i.e., load becomes higher). An interesting result that came out of our analysis is that the 
energy efficiency of the proposed scheme increases as the network becomes denser. Based on 
this result, we proposed the algorithm to select the relays that we referred to as interference- 
harnessing routing. In interference-harnessing routing, in contrast to the interference-aware 
routing, the routing criterion selects relays to exploit interference, instead to avoid it, thereby 
greatly increasing the energy efficiency in the network. By using interference-harnessing routing 
and the proposed QMF, the rate performance of multihop routing can be achieved by significantly 
reducing transmit power at each relay. 

Furthermore, to overcome the drawback of QMF/NNC which depends on prohibitively com¬ 
plex joint decoding at the destination, we developed a low-complexity stage-by-stage successive 
MIMO decoding. The proposed decoder can be implemented using a conventional MIMO 
decoding as ML, ZF, linear MMSE, and IF receivers. By comparing these various receiver 
architectures, we identified that, when the number of hops is less than 4, IF receiver provides 
a satisfactory performance within a I-bit from the ML performance and shows much better 
performance than multihop routing. In practice, the maximum number of hops is usually limited 
due to a delay constraint. For such case, the proposed scheme can provide a substantial gain over 
the multihop routing with similar encoding/decoding complexity, by implementing our scheme 
with IF linear receiver and a point-to-point channel code (e.g.. Turbo code, LDPC code, and 
Polar code). 
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